s-wave Cooper pair insulators and theory of correlated superconductors 
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The pseudogap state of cuprate high-temperature superconductors has been often viewed as either 
a yet unknown competing order or a precursor state to superconductivity. While awaiting the reso- 
lution of the pseudogap problem in cuprates, we demonstrate that local pairing fluctuations, vortex 
liquid dynamics and other precursor phenomena can emerge quite generally whenever fermionic 
excitations remain gapped across the superconducting transition, regardless of the gap origin. Our 
choice of a tractable model is a lattice band insulator with short-range attractive interactions be- 
tween fermions in the s-wave channel. An effective crossover between Bardeen-Cooper-Schrieffer 
(BCS) and Bose-Einstein condensate (BEG) regimes can be identified in any band insulator above 
two dimensions, while in two dimensions only the BEG regime exists. The superconducting tran- 
sition is "unconventional" (non-pair-breaking) in the BEG regime, identified by either the bosonic 
mean-field or XY universality class. The insulator adjacent to the superconductor in the BEG 
regime is a bosonic Mott insulator of Gooper pairs, which may be susceptible to charge density 
wave ordering. We construct a function of the many-body excitation spectrum whose non-analytic 
changes define a sharp distinction between band and Mott insulators. The corresponding "second 
order transition" can be observed out of equilibrium by driving a Gooper pair laser in the Mott in- 
sulator. We explicitly show that the gap for charged bosonic excitations lies below the threshold for 
Gooper pair breakup in any BEG regime, despite quantum fluctuations. Our discussion ends with 
a view of possible consequences for cuprates, where antinodal pair dynamics has certain features in 
common with our simple s-wave picture. 



I. INTRODUCTION 

The microscopic origin of high temperature supercon- 
ductivity, particularly in copper oxides, remains one of 
the greatest unresolved challenges in condensed matter 
physics. All systems known to superconduct at relatively 
high temperatures exhibit strong electronic correlations 
and great complexity, which the standard well-controlled 
theoretical techniques are ill-equipped to handle. Most 
novel approaches to correlated superconductors either fo- 
cus on spin dynamics^"— , or charge dynamics2rJ^ in their 
attempt to explain the multitude of unconventional prop- 
erties of the correlated normal states. This is naturally 
justified by proximity to either antiferromagnetic or su- 
perconducting phases. 

The charge-oriented approaches are based on the as- 
sumption that strong pairing correlations, or locally sta- 
ble Cooper pairs, can be held responsible for the pseudo- 
gap and normal state features. This "precursor" picture 
indeed finds support in a number of experimentsi^^— . 
If so, a prominent physical description of such normal 
states, originating in duality^", is a liquid of quantum 
vortices and antivortices, whose free motion dynamically 
disorders the phase of the complex order parameter and 
thus preemts "Bose condensation" of such locally stable 
Cooper pairs. 

The freedom in principle to either form or not form 
locally stable Cooper pairs opens a possibility of hav- 
ing two kinds of normal states. Conventional normal 
states have unpaired fermions, and are exemplified by 
Fermi liquids, band insulators and integer quantum Hall 
states (IQHS). The second kind are unconventional nor- 



mal states which involve manifestations of Cooper pair 
stability and exhibit strong correlations. Examples of 
such states discussed in literature are pair density waves 
(PDW)iiii^ and a class of fractional quantum Hall states 
(FQHS) with either an even-denominator filling factor or 
non-abelian statistics^i"— . 

In this paper we explore conditions for the develop- 
ment of pairing correlations in quantum insulators at 
zero temperature, and their observable consequences. 
Our analysis is focused on the nearly unitary scatter- 
ing in the s-wave channel, for several reasons. First, 
the theoretical convenience of studying the unitarity 
limit is well matched with the present-day cold atom 
experiments which routinely access unitarity in the s- 
wave channel. Second, demonstrating insulators dom- 
inated by s-wave Cooper pairs reveals that correlated 
non-superconducting states with pairing fluctuations are 
not the privilege of unconventional pairing mechanisms 
such as found in d-wave cuprates. An indirect impli- 
cation is that certain unconventional properties of nor- 
mal states in high- Tc superconductors could indeed be 
the consequence of certain universal aspects of pairing 
fiuctuations, rather than the specific microscopic mech- 
anisms responsible for pairing. In this manner, our re- 
sults indirectly support the "charge-oriented" or "vortex 
liquid" picture of cuprates, with an added benefit of tak- 
ing into account the fermionic excitations. The s-wave 
Cooper pair insulator may even be susceptible to charge 
density wave order as we shall briefly discuss, which is 
another interesting similarity to cuprates. Finally, by fo- 
cusing on the s-wave channel, we can sidestep the com- 
plications that are associated with gapless fermions in 
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nodal d-wave superconductors and the closely related is- 
sue of the absence of true "real space" pairs - and thus 
of the BEC counterpart to the BCS regime - in sys- 
tems where the pairing correlations are intrinsically re- 
pulsive. In those circumstances, strong quantum fluctu- 
ations ultimately lead to break-up of Cooper pairs and 
the low-energy physics of ensuing "d-wave duality" de- 
mands description which goes beyond the quantum liquid 
of vortex-antivortex excitations^!. Nevertheless, even in 
such systems, an intermediate state - having strong pair- 
ing correlations but lacking the long-range superconduct- 
ing order - is typically present in the phase diagram, and 
many of its properties will reflect the universal aspects 
of our s-wave insulator— 

A sharp statement which follows from renormalization 
group (RG)2S: and arguments in this paper is that precur- 
sor s-wave pairing in the insulating ground state is tied 
to the BEC regime which is separated from the weak- 
coupling BCS regime by a scattering resonance. Sim- 
ilar pairing phenomena have been studied in the deep 
quantum Hall limit^^ and at finite temperatures^L^^. 
Presently we view the BEC and BCS regimes as effective 
and appropriate for a particular band insulator or other 
gapped state of fermionic quasiparticles. Since the lat- 
tice potential frustrates the motion of microscopic parti- 
cles, the strength of microscopic interactions is effectively 
enhanced to pull the system toward the effective BEC 
limit. The RG indicates that a BEC-BCS crossover with 
attractive interactions can be observed only above two 
dimensions, where the BEC regime requires a sufficiently 
strong coupling. In two dimensions, however, even ar- 
bitrarily weak coupling puts the system into the BEC 
regime, implying the existence of charged bosonic exci- 
tations in the vicinity of the superconducting transition. 
The second order superconductor-insulator transition in 
a BEC regime is captured by a purely bosonic effective 
theory, as the fermionic excitations remain gapped across 
the transition. Therefore, the insulating state is a bosonic 
Mott insulator, which can be also viewed as a relativistic 
quantum liquid of vortices and antivortices if the tran- 
sition is driven at a fixed quasiparticle density in the 
grand-canonical ensemble. 

A clear and observable manifestation of a bosonic Mott 
insulator adjacent to the superconducting phase is the 
bosonic mean-field or XY universality class of the tran- 
sition. However, a Mott insulator should then have a 
more robust identity, something to distinguish it from a 
band insulator which can be turned into a superconduc- 
tor only by closing the fermion gap. Here we construct 
a special spectral function which captures non-analytic 
changes of the many-body spectrum associated with the 
disappearance of excitations with bosonic statistics. This 
function behaves as an order parameter at a second order 
transition when the insulator changes its character from 
Mott to band. However, no thermodynamic second or- 
der phase transition takes place, since only the excitation 
spectrum evolves in a non-trivial way. This kind of sharp 
changes cannot be observed in equilibrium at zero tem- 



perature, but could appear as actual phase transitions if 
the system is driven out of equilibrium. 

A concrete non-equilibrium distinction between 
bosonic Mott and band insulators that we propose is the 
ability to function as a Cooper pair laser. A coherent ex- 
cited state which carries a macroscopic charge or current 
is possible to generate and sustain in a Mott insulator, 
but not possible in a band insulator. The spectral func- 
tion we construct is a representation-independent order 
parameter of the Cooper pair laser, as this driven con- 
densate spontaneously breaks the "U(l) symmetry". 

An important part of our analysis is a proof that in- 
deed a two- or three-dimensional insulator in an identifi- 
able BEC regime is characterized by the existence of (i) 
low energy bosonic excitations and (ii) gapped high en- 
ergy fermionic excitations. Specifically, the fermion gap 
remains finite across the transition at which the boson 
gap vanishes. We show this using a model band insu- 
lator of quasiparticles and quasiholes near unitarity as 
a starting point. We consider both the non-relativistic 
and relativistic effective bosonic dynamics. The propa- 
gator for bosonic excitations is calculated in a perturba- 
tive random phase approximation manner from the ex- 
act fermion propagators and effective short-range inter- 
actions embodied in a renormalized vertex function. The 
mathematical forms of these propagators and vertices are 
qualitatively fixed by the fact that the ground state is 
an insulator. In this sense, we derive exact conclusions 
about the spectra of bosonic and fermionic excitations, 
taking into account all fluctuation effects. These findings 
complement the recent RG analysis^^ which was not able 
to access states with locally stable Cooper pairs. 

While cuprates are fundamentally different than the 
s-wave model we analyze, certain remarkable similarities 
between them are very much worth exploring. If one fo- 
cuses on the antinodal regions in the Brillouin zone, the 
superconducting transition in underdoped cuprates looks 
just like the superconductor-insulator transition we dis- 
cuss in the two-dimensional BEC regime. The pairing 
gap is indeed rooted in the antinodal regions where short- 
range spin correlations alone would gap-out fermionic ex- 
citations. The existence of gapless nodal quasiparticles, 
which prohibits pair bound states and thus spoils the full 
BEC analogy, is almost a secondary effect of the d-wave 
pairing. What is not spoiled is the dominance of bosonic 
fluctuations in the pseudogap state, owing to the linearly 
vanishing density of states of the nodal quasiparticles. 

This paper is organized as follows. In section we 
define the model of a band insulator with attractive in- 
teractions and discuss its insulating states. Following the 
derivation of low energy effective theories for the insulat- 
ing states of our model, we explain in detail the sharp dis- 
tinction between Mott and band insulators and ways to 
observe it out of equilibrium. In section|nT]we present an 
estimate of the Cooper pair insulator regime in the phase 
diagram of a realistic two-dimensional cold-atom system 
using a numerical mean-field calculation. The section 
IIVI deals with the stability of Cooper pair insulators and 
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presents the perturbative fluctuation analysis. Conclu- 
sions are summarized in the discussion section, including 
experimental prospects and a more detailed qualitative 
view of the pseudogap physics in cuprates. 



II. BOSE V.S. FERMI INSULATOR 

The most promising circumstances for the occurrence 
of strongly correlated states can be found whenever there 
is a macroscopic degeneracy of states which can be lifted 
even by weak perturbations. One can ask what happens 
to an integer quantum Hall insulator when electrons ac- 
quire attractive interactions. Sufficiently strong attrac- 
tive interactions in this topological band insulator can 
yield Cooper pairing and re-entrant superconductivity 
at zero temperature. However, direct BCS-like pairing 
(or pair-breaking) transitions are preempted by first or- 
der transitions and masked by intervening strongly cor- 
related insulating phases^^. Such correlated insulators 
are obtained from superconductors by quantum vortex 
lattice melting, and hence can be called "vortex liquids" . 
While their properties are generally not universal, their 
existence at zero temperature provides a sharp thermo- 
dynamic distinction between paired and unpaired non- 
superfluid states when time-reversal symmetry is vio- 
lated. 

Here we explore whether some sharp distinction be- 
tween paired and unpaired insulators can be found in 
generic time-reversal invariant circumstances. A band 
insulator with attractive interactions is indeed suscep- 
tible to symmetry breaking involving phase-incoherent 
Cooper pairs. A recent study has shown that pairing in- 
stabilities in ideal multi-band insulators always occur at 
finite rather than zero pair-momenta. Such instabilities 
produce supersolid pair density wave (PDW) states^S,. 
The pairing wavevector is generally incommensurate with 
the lattice as a consequence of the non-trivial crystal- 
momentum dependence of vertex functions for attractive 
contact interactions. While fluctuations cannot elimi- 
nate the finite-momentum instability, a simple argument 
shows that they could have a profound influence on se- 
lecting the wavevector. The incommensurate supersolid 
PDW, which naively results from this instability, would 
be highly frustrated in the lattice potential. The result- 
ing vast multitude of degenerate PDW ground-states can 
be easily mixed by quantum fluctuations until the U(l) 
symmetry is restored, or a supersolid at a commensurate 
wavevector is stabilized. We will show in this paper that 
strong quantum fluctuations near the second order tran- 
sition can result with a bosonic Mott insulator of Cooper 
pairs which intervenes between the superconductor and 
band insulator states in the phase diagram. Then, the 
fluctuations which stabilize a commensurate order in the 
supersolid phase are expected to do the same in the Mott 
insulator, and produce a commensurate insulating PDW. 
The obtained insulating phase is thermodynamically dis- 
tinct from any featureless band insulator due to trans- 



lation symmetry breaking. It is also distinguished from 
ordinary "single-fermion" charge density waves by its low 
energy charged bosonic excitations. Translational sym- 
metry breaking is possible even though the number of 
Cooper pairs per site is an integer, due to the orbital 
degrees of freedom carried by Cooper pairs, and peculiar 
non-analytic features of bosonic propagators and vertices 
at the zero momentum^^. 

Most generally, however, the insulating state proxi- 
mate to the superconductor need not have broken sym- 
metries, or any other thermodynamic distinction from 
the ordinary band insulator. Nevertheless, it is still pos- 
sible to sharply distinguish it from the band insulator 
in out-of-equilibrium situations. We will argue that our 
model indeed supports an insulator with infinite-lifetime 
Cooper pairs. The physical picture is as follows. While 
the fermionic excitations have a large bandgap, bosonic 
collective excitations must have a much smaller gap in the 
insulator near the second-order superconducting transi- 
tion. The conservation of energy, momentum and charge 
protects the inflnite life-time of these bosons (in an iso- 
lated system at T = 0) despite interactions among them. 
While some relaxation processes can take place result- 
ing with decays into fermion pairs, excited momentum 
and charge carrying collective modes below the bandgap 
cannot completely vanish without coupling to the envi- 
ronment. Therefore, by coupling to a suitable external 
driving system, it is possible, as a matter of principle, to 
excite this bosonic insulator into a macroscopically co- 
herent state which spontaneously breaks U(l) symmetry. 
Such a driven (non-relativistic) condensate is analogous 
to the state of (relativistic) photons produced by a laser. 
If the boson gap is pushed above twice the fermion gap, 
for example by reducing the strength of pairing inter- 
actions, then the coherent bosonic excitations cease to 
exist and the boson-laser non-equilibrium state becomes 
impossible. 

Formally, we deflne phase transitions as points in the 
parameter space {fj,i,i>i} at which the thermodynamic 
potential J-{fii^ Vi) is not an analytic function. If we know 
the Hamiltonian H{fii), we can construct the density ma- 
trix p as a function of both microscopic parameters fii and 
external parameters i>i which can describe the environ- 
ment (temperature, details of the non-equilibrium drive, 
etc.). In thermal equilibrium at temperature T = 1/(3 
we simply have p = z^^e"^^ , where z — ti{e~^^), and 

— —Tlogz. Clearly, at T = only the ground-state 
manifold of the Hamiltonian enters J-, so that quan- 
tum phase transitions are possible only if the ground- 
state spectrum undergoes non-analytic changes, such as 
changes of degeneracy. However, the density matrix at 
finite temperatures, or in generic non-equilibrium situa- 
tions (even at T = 0) can inherit non-analytic behavior 
from the excited state spectrum. Deriving thermody- 
namic functions in these general cases provides mathe- 
matical routes for defining "non-equilibrium phase tran- 
sitions" and sharp distinctions between states which may 
be thermodynamically equivalent. The Cooper pair in- 
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sulator described above is an example of this principle. 



A. The model and insulating states 

As a generic model of a band insulator let us consider 
the imaginary-time action of neutral fermionic particles 
with attractive interactions U, in a lattice potential V{r): 



S= dT 



dVd''r2f/(|ri - r2|)V'j(ri)V't(ri)^I(r2)V'^(r2) 



This is a microscopic multi-band model defined in contin- 
uum space and not a priori tied to the vicinity of a critical 
point. Let us imagine that we can control the chem- 
ical potential /i, lattice V{r) and the overall strength 
of attractive interactions U between fermions. In cold- 
atom realizations of the model it is easy to manipulate U 
through Feshbach resonances and the depth of an opti- 
cally generated lattice potential, while condensed matter 
systems allow direct control of the chemical potential. 

Suppose we place the chemical potential somewhere in 
a bandgap of the fermion spectrum and start with strong 
attractive forces for which the ground state is a super- 
fluid. Then, we gradually reduce U at zero temperature. 
By focusing on three dimensions and not too shallow 
lattice potentials we can ensure the following changes. 
First, the system undergoes a transition to an insulator 
at some critical U = Ugf ■ Goldstone modes of the super- 
fluid acquire a gap, but initially this gap is much smaller 
than the gap for fermionic excitations. In this regime, 
the low-energy effective theory can be obtained from ([T]) 
by applying a Hubbard- Stratonovich transformation in 
the particle-particle channel to decouple the interaction, 
followed by the fermion field integration. The lattice the- 
ory which describes the superconducting transition, de- 
rived in Appendix|3 contains only charged bosonic fields 
(Cooper pairs): 



(2) 
off 



dr 



d_ 



ijkl 



(2) 

The indices label both the lattice sites and orbital 

degrees of freedom, and -fCj'j^[£] is an operator function of 
X with a well defined Taylor expansion about x ~ 0. Note 
that all terms in this action have local character owing 
to the finite gap of the fermionic fields whose fiuctua- 
tions produce the effective boson dynamics. This means 
that the momentum p and frequency dependence of 
all couplings is analytic at p = and = 0, so that the 
bosonic fields b correspond to physical bosonic excitations 
even though their mathematical origin is in a Hubbard- 
Stratonovich transformation. The insulator adjacent to 
the superconductor is essentially a bosonic Mott insula- 
tor. 



As U is decreased even further, the mass terms for all 
boson fields grow. The boson gap increases and the en- 
ergy range spanned by the long-lived bosonic excitations 
keeps shrinking, being limited from above by the thresh- 
old for decay into fermions. At some other value U = Um 
the boson gap becomes equal to twice the fermion gap. 
Beyond that point, the lowest energy excitations are not 
bosonic, and the above effective theory is not applicable. 
Instead, the theory which captures dynamics at the low- 
est energies can be obtained by writing the multi-band 
lattice version of the single-channel model ([T]) , or formally 
by integrating out the b fields in (jA9p or (|A12[) : 



^^ff ^ - / dr 



E 



d_ 



.ft f 



^^3klfi^fj-\fkifli 

ijkl 



(3) 



This theory describes a band insulator. 

The theory ([2]) has only bosonic excitations at low 
energies, while ([3]) has only fermionic excitations. The 
change of statistics implies that the many-body spec- 
trum cannot smoothly evolve between these two limits. 
One way to see the non-analytic change is from the spec- 
tral function of "charged" bosonic excitations, analyzed 
in section HVl which loses a sharp coherent peak between 
the two limits. However, we shall pursue a different route 
in the next section, and consider the density of many- 
body states in order to discuss implications for out-of- 
equilibrium dynamics. 



B. A sharp distinction between Mott and band 
insulators at T = 

Let us now start from a bosonic Mott insulator in the 
limit < Ab ^ 2At, where Ab and Af are the gaps 
for bosonic and fermionic excitations respectively. The 
elementary low energy excitations are gapped bosonic 
modes labeled by crystal wavevector k and an orbital in- 
dex which will be suppressed. The many body spectrum 
is further characterized by integer numbers of quanta 
Nk > excited in different bosonic modes. 

Since all effective interactions between bosons are me- 
diated by gapped fields, they are short-ranged in real 
space. Consequently, interaction effects become notice- 
able in the many-body spectrum only when the density of 
bosonic excitations becomes finite, that is when the aver- 
age real-space distance between them becomes finite. As 
long as Nk is not macroscopically large we can write 
the following expressions for the total energy E and mo- 
mentum P: 



E 



(4) 



The individual mode energies are > Ab, while the 
mode momenta p^. can be arbitrarily small. The total en- 
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ergy is conserved, and the total momentuni is conserved 
modulo first Brillouin zone. 

Consider the many-body density of states (DOS) as a 
function of conserved quantities: 

p(£;,p) = i^5(£;-i?„)<5(P-p„) , (5) 

n 

where n labels all many-body Hamiltonian eigenstates, 
and V is the d-dimensional system volume. Let us define 
a modified "single-mode" DOS: 

p'(k) = hm hm — / dSe [ d'^Sp 

A, A^ 

N 

The small energy and momentum-space volumes, and 
Ap respectively, should be taken to zero only after taking 
the thermodynamic limit V — > oo. This ensures that the 
above energy and momentum integrals of the many-body 
DOS (O converge in the thermodynamic limit, as the 
number of modes in a fixed momentum-space volume Ap 
scales as V. The sum over N extracts the essential prop- 
erty of the many-body spectrum, the extent over which 
the spectrum has the "period" e^. The single- mode DOS 
p' would have been divergent even in a finite volume in 
the absence of interactions between bosons. However, 
as already pointed out, interaction corrections to spectra 
scale as 1/ V at a fixed density (or N) , so that finite devia- 
tions of the many-body spectra from the non-interacting 
form ^ a.t N (X V cut off the divergence of p' . The 
factor of 1/V in © then keeps p' finite in the thermo- 
dynamic limit. Therefore, p' will by construction reflect 
the density of bosonic quanta which can be excited in a 
single mode before interactions become too strong, which 
suggests an "order parameter" interpretation. We shall 
immediately see that the behavior of p' indeed resembles 
a second order transition at the crossover between Mott 
and band insulators. 

Tuning the Hamiltonian parameters to increase the 
bosonic gap results in the gradual reduction of p' . The 
excitation density needed to visibly modify the many- 
body spectrum away from the free-particle form becomes 
smaller, as the spatial range of effective interactions be- 
tween low-energy bosons grows. This can be seen from 
the Feynman diagram in Figlljb) which captures the two- 
body interactions between bosonic excitations. Let us 
focus on low energy bosons in the Mott insulator whose 
decay into fermions in a collision is prohibited by momen- 
tum and energy conservation. Evaluating this diagram 
(see appendix |B]) at zero energy transfer (po = 0) and 
finite momentum transfer (p = |p| ^ 0) in the limit of 
low energy and momenta of incoming particles reveals 
a scale = 4m(2Af — Ab) which appears alongside 
p^. The diagram is convergent in the p — > limit when 
r is finite, but becomes infra-red divergent for p — > 
when Ab = 2Af. Therefore, r acts as a characteristic 




FIG. 1. Diagrammatic representation of (a) quadratic and 
(b) quartic terms in the effective action ((2]) describing the 
dynamics of bosonic excitations. Solid lines represent fermion 
propagators G(fc^) — {uj—£k)~^, while the faded external lines 
represent Cooper pair fields. The diagram (a) generates the 
inverse pairing susceptibility n(g^) given by p2|) and (|17|) . 
while (b) generates effective interactions between two bosonic 
excitations. As long as the boson gap Ab is smaller than twice 
the fermion gap Af (implicit in Sk), the poles of low energy 
bosons obtained from n(g^) — occur at real frequencies, and 
the interaction (b) between low energy bosonic excitations 
has a finite range r ~ [4m(2Af- Ab)]"^/^ Beyond Ab > 2Af 
the interactions would have been effectively long-ranged if the 
boson poles survived. However, the boson pole is lost in such 
a BCS regime (see section HV)) . 



length-scale over which the short-range effective interac- 
tions between bosonic excitations are extended. 

When the boson gap reaches and exceeds twice the 
gap for fermionic quasiparticles, the sum over iV in ([S]) 
gets cut-off at an N/V — > 0, so that p' — > in the 
thermodynamic limit for every single-particle mode. If 
bosonic excitations could survive, the inter-boson inter- 
actions would become effectively long-ranged (r — )■ oo) 
and would completely remove the period from the 
many-body spectrum. However, more importantly, all 
bosons, including the lowest-energy ones, now inevitably 
decay into fermion pairs regardless of whether they col- 
lide with other bosons or not. We shall find in the 
section IIVI that bosonic poles must completely disap- 
pear when Ab > 2Af, resulting in a BCS regime. The 
only elementary excitations are fermionic quasiparticles 
whose spectrum cannot admit multiply occupied single- 
particle states. The insulator in this regime is a band 
insulatori^Si^. 

Therefore, p'{k — 0) is an analytic function of the 
many-body spectrum in a finite-size system, but under- 
goes a non-analytic change in the thermodynamic limit 
at the crossover Ab = 2Af between a bosonic Mott 
and fermionic band insulator. This non-analytic feature 
stems from the fundamental change in the many-body 
spectrum of excited states, and can be observed as a 
"second order phase transition" in out-of-equilibrium cir- 
cumstances. 

A similar non-analytic change occurs in p{k) for any 
fc when Cfc = 2Af. The abundance of singularities at 
different k does not trivialize the matter since the en- 
semble averaged operators can pick up different sets of 
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these singularities in different out-of-equilibrium circum- 
stances. For example, p{k) can be regarded as an or- 
der parameter for a driven condensate at wavevector fc, 
which we discuss in the following subsections. The quan- 
tity p{k = 0), however, is special because its non-analytic 
behavior is the last to disappear when the boson gap is 
gradually increased, in addition to fc = representing the 
most symmetric excited state. At the equilibrium phase 
transition to the superfluid (Ab — >■ 0), p{k = 0) becomes 
proportional to the standard order parameter associated 
with U(l) "symmetry breaking" . 



C. Cooper pair laser 

A finite p'(fe) in ([6]) is an indication of coherent many- 
body excitations. Here we argue on more physical 
grounds that the existence of such states which carry 
macroscopic charge or current is indeed protected by con- 
servation laws and energetics in the insulating phase with 
short-range interactions. Therefore, at least in principle, 
a bosonic Mott insulator can be turned into a Cooper pair 
laser by some appropriate external drive. The coherent 
quality of bosonic low energy dynamics is experimentally 
observable at least out of equilibrium. 

Consider preparing an ensemble of A'^b coherent Cooper 
pairs, all occupying the same excitation mode with mo- 
mentum p and energy ep. Let there be no external par- 
ticle reservoir or heat bath in contact with the system. 
Due to interactions this is not an eigenstate of the many- 
body spectrum, so particles will be scattered out of the 
coherent wave. Some of the Cooper pairs can gain enough 
energy to break into fermion pairs. After many scatter- 
ing events, a steady state contains bosons and Nf 
fermions, where a finite fraction of the remaining bosons 
may still participate in a dressed coherent wave. 

Being in a bosonic Mott insulating state with Ab < 
2Af ensures the existence of infinite-lifetime single- 
particle bosonic excitations with energy Ab < e < 2Af. 
However, two such low-energy bosons can still have 
enough total energy to leave behind a pair of free fermions 
after a collision as one of the bosons is pushed to high en- 
ergy above 2Af and the other drops to a lower energy and 
takes a recoil momentum. Hence, we must first explore 
whether at least some bosons from the initial ensemble 
can survive collisions. The conservation of energy and 
charge requires: 



A^bCp 
2Ab 



(7) 



2Ni 



Here > and Kf > are the amounts of kinetic 
energy per surviving low-energy boson and excited high- 
energy fermion respectively, some of which can be con- 
sidered heat. It follows that 



K 



£p -2{Af + Kf) 
Ab-f ifb-2(Af + iff) 



Since N^^ < N\^, we must have Cp > Ab + K\^- If the 
energy ep of the initial Cooper pairs is larger than the 
threshold 2Af for the decay into two fermions, then Ki 
can grow until becomes zero, implying the decay of all 
Cooper pairs. The increase of entropy likely leads to this 
outcome. Otherwise, ep < 2Af makes the numerator in 
([8|) negative (for any ATf > 0), so that the denominator 
must be negative as well in order to keep A^^ positive. 
Then, a consequence of ATb > is 



K 



> 



1 - 



2(Af+j<-f) 
2(Af+A:f) 



(9) 



(8) 



In other words, A'^ is a finite fraction of Ab, unless ep 
reaches or exceeds 2(At + Kf) when A'^ can drop to zero. 
Therefore, if a bosonic mode with energy Ab < ep < 2Af 
is macroscopically excited, then the resulting steady state 
will still have a macroscopic number of surviving bosons. 

The conservation laws merely protect the ensemble 
against losing all bosons to fermions. However, we must 
yet argue that a finite fraction of the remaining bosons 
participates in a macroscopic coherent wave. This is in 
fact straight-forward and on the same footing as the case 
for the stability of equilibrium condensates. Let us boost 
to the "center-of-mass" reference frame of the ensemble 
where a coherent wave would look like a regular static 
condensate. If no bosons were to remain in the coherent 
wave due to collisions, then they would have to be local- 
ized into a "crystal" in this frame at distances a oc n~d 
from each other, where n is their density. By Heisen- 
berg uncertainty principle the localization (kinetic) en- 
ergy per particle would scale as a~^, while the potential 
energy due to short-range interactions can be modeled 
as an exponential-like function e"*"*. Clearly, regardless 
of the numerical co-cfficicnts in these functional forms, 
localization becomes too costly at large enough a, so the 
"crystal" melts at a finite low density n. The neglect of 
the lattice potential for the purposes of this argument is 
justified in the limit of n much smaller than one particle 
per site. 

Therefore, coherent many-body excitations are stable 
if they hold a macroscopic but not too large number of 
Cooper pairs. Relaxing the ensemble to a maximum en- 
tropy state with the given macroscopic momentum (by 
the external Cooper pair laser pump) would leave behind 
a macroscopic number of bosons in a coherent wave be- 
cause this minimizes the energy of particles which do not 
take away heat. 

Another possible issue are finite wavevector dynamical 
instabilities of the excited coherent wave^^r^^^. However, 
we do not expect such instabilities to arise near conven- 
tional second order superconductor-insulator transitions 
that we focus on. Since the equilibrium Mott insulator 
has an integer number of Cooper pairs per site, its lowest 
energy excitations live at the ordering wavevector(s) of 
the adjacent superconducting state (uniform or PDW). 
Any roton excitations at other crystal wavevectors are 
separated by an additional gap, so they do not provide a 
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channel for dynamical instabilities at slow excited flows. 
Note that the Cooper pair laser is fundamentally dynam- 
ically unstable at zero^'^'^ or any PDW wavevectors - 
this is precisely what the "laser pump" has to defeat in 
order to drive a coherent condensate. 



D. A non-equilibrium T — phase diagram from a 
hypothetical experiment 

While many out-of-equilibrium phase transitions have 
been considered in literature, the overall theory is far 
less developed than in the case of equilibrium (thermo- 
dynamic) phase transitions. Systematic studies of criti- 
cal dynamics, expressed through properties such as re- 
laxation rates and transport, elucidated a number of 
non-equilibrium universality classes^. In many cases 
critical dynamics was analyzed in the vicinity of equi- 
librium critical pointsS^Ti^. Also, very detailed stud- 
ies of non-equilibrium physics have been devoted to 
one-dimensional systems^, where an additional incen- 
tive comes from integrable cases which lack ergodicity 
and do not reach thermal equilibrium if decoupled from 
the environment. This interest in one-dimensional sys- 
tems has recently grown in the context of ultra-cold 
atoms^^— . On the experimental side, various demon- 
strations of phase transitions caused by driving a system 
out of equilibrium have been demonstrated^^—. A more 
complete phase diagram is parametrized by both equilib- 
rium and non-equilibrium parameters, so that the same 
transition can be driven either in equilibrium or out of 
equilibrium-- . 

Here we attempt to put the previous discussion in the 
context of a fictitious experiment which aims to probe the 
boson-dominated dynamics in the insulating state. We 
shall imagine that a laser light of wavelength A which can 
excite collective bosonic modes above the ground state is 
shined on the sample. For this purpose we are interested 
in absorptions of an even number of photons in which 
there is no angular momentum transfer (the laser should 
be linearly polarized) . As long as the momentum of ab- 
sorbed photons is negligible, excitations beyond the ab- 
sorption edge must contain more than one bosonic mode. 
Single photons can also be absorbed, but the threshold 
for such processes is given by twice the fermion gap and 
lies above the threshold for the two-photon absorption 
in the interesting Bose insulator regime. The schematic 
phase diagram is shown in Fig (21 

If the ground-state is a superfluid (SF), any laser wave- 
length A will result in some partial depletion of the 
ground-state condensate and promotion of some Cooper 
pairs to finite energies. If the energy of promoted Cooper 
pairs is larger than twice the fermion gap 2Af , the excited 
pairs decay and leave behind a population of free fermions 
which can conduct current as a metal. Otherwise, the 
promoted Cooper pairs can survive as long-lived low en- 
ergy excitations and form a coherent state, a Cooper pair 
laser (CL). 
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FIG. 2. A schematic zero-temperature phase diagram of the 
model (m driven out of equilibrium by a laser of wavelength 
\ (c = h = l). We assume that the laser induces at least 
two-photon processes in which Cooper pairs are promoted 
to higher energy states. In equilibrium at = there 
are two thermodynamic phases, superfluid (SF) and insulator 
separated by a critical point (filled circle). The many-body 
spectrum undergoes another non-analytic change in the ther- 
modynamic limit (open circle) which affects excitations and 
introduces a sharp distinction between a bosonic Mott insu- 
lator (MI) and a fermionic band insulator (BI). This "tran- 
sition" may not be thermodynamic, there is no latent heat, 
onset of symmetry breaking or topological order. Going out 
of equilibrium allows dynamic population of higher energy 
states at expense of depleting the ground state. If the ex- 
cited states have bosonic character, light absorption can yield 
a driven condensate, or a Cooper pair laser (CL). However, if 
the two-photon energy exceeds twice the fermion gap, the ex- 
cited Cooper pairs decay into free quasiparticles and a driven 
metal (M) is obtained. When the ground-state is a superfluid, 
the driven excited state can coexist with the depleted ground- 
state superfluid (SF-I-CL and SF-I-M). The labels Ab and Af 
denote gaps for bosonic and fermionic excitations, which are 
controlled by the chemical potential ^. The vertical dashed 
line indicates the non-analytic change of the spectrum as a 
function of only the equilibrium parameters, which may not 
be felt within the ground-state manifold. 



If the ground-state is an insulator, it may still be pos- 
sible to obtain a Cooper pair laser, as long as the excited 
low-energy Cooper pairs cannot decay into fermions. 
This possibility exists if the boson gap Ab is smaller 
than twice the fermion gap 2Af. Otherwise, photon 
absorption yields an excited population of interacting 
but free fermions. The insulating ground-state which 
can be excited into a driven condensate is designated 
a bosonic Mott insulator (MI) in Figl^l The analy- 
sis in the previous subsection points to the conclusion 
that the many-body spectrum of the equilibrium Hamil- 
tonian undergoes a qualitatively significant non-analytic 
change when Ab — 2Af, the open circle in Figl5] The 
spectrum remains gapped and it is not clear whether 
the non-degenerate ground-state energy experiences any 
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non-analytic dependence on the equilibrium parameters. 
Therefore, we cannot claim any thermodynamic phase 
transition between the Mott and band insulators, even 
though formally there exists a sharp distinction between 
them, visible in the structure of excited many-body 
states. 

Hence, we argue that the phase diagram contains a 
special MI-BI crossover, depicted by the dashed vertical 
line in Fig 121 which has some sharp signature in the spec- 
trum and could be physically described as pair-breaking. 
The envisioned hypothetical experiment could only indi- 
rectly sense the existence of this crossover, by probing 
the possibility of creating a Cooper pair laser from the 
insulating equilibrium state. 

Realizing Cooper pair laser states in realistic systems is 
of course very challenging. While dissipation due to ab- 
sorption can always be defeated by an appropriate laser 
pump which takes advantage of the stimulated emission, 
it remains unclear what could be an appropriate pump 
in a given system, and whether it is feasible in practice 
(probably not in solids). In trapped cold atom gases, 
it is at least possible to make an ordinary "condensate" 
of bosons in an optical lattice, then suddenly change the 
lattice depth to create conditions for an equilibrium Mott 
insulator at the trap center. This system, initially in a 
"laser" state, would take a long time to relax into the 
many-body ground state. We just do not know how to 
continuously pump this atomic Cooper pair laser. 



The bubble diagram in the particle-particle channel 
from FigUfa) 



n 
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is the first step toward finding a superfluid instability. 
Here k and q are conserved first Brillouin zone (BZ) 
wavevectors, n and m are "band-index" quantum num- 
bers for bosonic and fermionic excitations respectively, 
£mk are bare fermion energies appropriate for the Bloch 
states (relative to the chemical potential), /(e) is the 
Fermi-Dirac distribution function, and F are vertex func- 
tions which can be derived from the fermionic ipmk and 
bosonic Bloch-state wavefunctions: 



C%,;,.2fc2-= / d^^1';qMV'.„,fc,(r)VW2W • (13) 



Since crystal momentum is conserved, linq-n' q' {^) = 
^nn'iQ, ^) X (27r)''^(<7 — q'). The expression (|T^ is ultra- 
violet divergent and we must regularize it by removing 
the contribution of vacuum zero-point fluctuations occur- 
ring at arbitrarily short length-scales. After regulariza- 
tion we may write the physical propagator Dnn' (q, ^) 
for the charged bosonic excitations: 



III. NUMERICAL RESULTS FOR A MODEL IN 
THE UNITARITY REGIME 



H„„'(q,l^)-H 



(14) 



Here we estimate the conditions for a bosonic insula- 
tor in a specific model, using a mean-field approxima- 
tion. We consider a model of neutral fermionic ultra- 
cold atoms in a lattice potential, whose interactions are 
tuned close to the Feshbach resonance. The imaginary- 
time "two-channel" action describing this system is: 



drd'^r 



47r 



l$l 



<l>^|Vj + h.c. 



(10) 



Only contact interactions are relevant near the unitar- 
ity fixed point, and they are decoupled by the complex 
Hubbard-Stratonovich field $. We shall use the simplest 
square lattice potential 



V{r)=2Vj2' 



2T:Xi 
aL 



(11) 



where Xi G {x,y,z...}, is lattice spacing and V is 
a tunable lattice amplitude. Integrating out the fermion 
fields ipa yields an effective action for the charged bosonic 
excitations represented by <&. 



where H^eg is the smallest eigenvalue of H„„/(qr, Q) eval- 
uated in vacuum at T = and q = 0, fi = 0. The 
formal calculation of H„„'(q, fi) in this expression must 
involve an implicit momentum cut-off A for internal mo- 
menta ki, but after regularization (|14|) is convergent in 
the A — cxD limit. The parameter v can be related to the 
scattering length a in two-fermion collisions in vacuum; 
in three dimensions, v = —1/a. 

The smallest eigenvalue n{fl, q) of ~D^^, (il, q) reveals 
superfluid instability when it becomes negative at zero 
frequency. Since 7r(f2, q) is independent of the representa- 
tion for Z?, we can choose any complete set of boson wave- 
functions ^riqi'i') which treats q as a crystal wavevector 
in the lattice (HH). In other words, ^nqi'i') can be Bloch 
states of any periodic potential with the same symmetries 
and lattice spacing as (|lip . and it will be convenient to 
carry out calculations using the plane- wave Bloch states, 
as if 1^ = 0, where n = G (a reciprocal lattice vector) 
and $Gq(r-) = e''^'^+"^'' . 

In the insulating phase 7r(0, q) > 0, but the boson gap 
may be small. We estimate the boson gap At from the 



fact that 7r(f2, q) « tt, 



(0) 



(2)n2 



51 for small boson 



energy f2, by noting that 7r(ri, q) — yields the dispersion 
relation for bosonic excitations. In fact, there are two 
bosonic modes, particle-like with gap Abe and hole-like 



9 



0.30 
0.28 
0.26 
0.24 
0.22 
0.20 
0.18 
0.16 




V 




(a) 



(b) 




FIG. 3. (color online) The phase diagram of neutral fermions with near unitary interactions, in a two dimensional lattice 
potential at T = 0. Going from weak to strong coupling: (a) u — 0.5, (b) f = 0, (c) u = —0.5 {Eq = fi^ /2ma\). The red thick 
line shows the boundary between the superfluid (SF) and insulating (I) phase. The thin red line outlines the first bandgap 
of bare fermions, toward which the SF-I transition line converges in the vanishing coupling limit v — ^ +oo. The shaded blue 
regions on the insulating side near the transition represent the difference between the boson and twice the fermion gap, darker 
colors corresponding to larger differences. These regions which grow with the strength of pairing are Cooper pair insulators. 
The dashed blue lines in subfigure (a) are trajectories in the parameter space along which transitions dominated by particles 
(p), holes (h), or both (ph) can occur. 



with gap Abv: 



persist beyond the mean-field approximation. 
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The insulator is bosonic if Abe < 2Ac or Abv < 2Av, 
where Ac and Av are fermionic particle and hole gaps 
which we read out from the band-structure at any 
given chemical potential in this mean-field approxima- 
tion. This estimate breaks down in the superfluid state. 

Figure 13] shows the phase diagram of a two-dimensional 
system in the vicinity of the superfluid-insulator tran- 
sition with the average density of two atoms per lat- 
tice site. The mean-field approximation finds a band- 
insulating state sitting inside the first bandgap when the 
pairing interaction is not strong enough. However, close 
to the superfluid boundary there is a region in which 
the gap for collective bosonic excitations is smaller than 
twice the fermion gap. The thickness of this Cooper pair 
insulator is small in the weak coupling limit, but grows 
with the strength of pairing. The spike feature in the 
shape of the Cooper pair insulator formally reflects the 
non-analytic dependence of the fermion gap on chemical 
potential (fermion gap is the distance between the chemi- 
cal potential and the nearest band-edge for a given lattice 
depth V). Physically, both particles and holes shape the 
bosonic dynamics in the vicinity of spikes, while only par- 
ticles or only holes are relevant near the band-edges. We 
cannot guarantee that such pronounced features would 



IV. THE STABILITY OF THE 
BOSE-INSULATOR 

In this section we explore the conditions for having 
gapped fermionic excitations across the T = superfluid 
transition. We are particularly interested in the possi- 
bility of "unconventional" superfluid transitions across 
which fermionic excitations remain gapped. The conven- 
tional BCS pair-breaking transition has gapless fermions 
in the normal state, at least just after the transition. 

We adopt a simple model of a band insulator in which 
there are two very broad bands, conduction and valence, 
separated by a direct bandgap Eg. Generalizations to 
other bandwidths and types of bandgap are straight- 
forward and are not expected to alter the main conclu- 
sions. The chemical potential residing in the bandgap 
controls the particle and hole gaps Ac and Av respec- 



tively, under the constraint Ac 



Eg. The particle 



energies in the conduction (c) and valence (v) bands are 
given by 



Sc(fc) = Ac 



2mc 



E^{k) = -Av - 



2mv 



(16) 



It is convenient to carry out calculations in the two- 
channel model pH)) . Restricting ourselves to zero tem- 
perature, the task of calculating the bosonic propagator 
D from ([14]) in the present model is not complicated. 
The numerator in (jl2p is non-zero at T = only if the 
two fermions in the bubble diagram are both above or 
both below the chemical potential. This presently means 
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that only mi = m2 = c and mi = m2 = v contribute. 
Assuming a very short-range microscopic interaction po- 
tential between fermions, we can ignore the momentum 
dependence of the vertex functions and write them as 
^mim2- Finally, we can choose the representation for 
IIiiii' which diagonalizes it and focus only on the small- 
est eigenvalue which captures the onset of superfluidity. 
There is only one set of n corresponding to the softest 
bosonic mode which we ought to focus on, leaving us 
with only two relevant dimensionless values that the ver- 
tex function can take, 7c = T^c and 7v = Fw- 

The Green's function D{q^) for the softest bosonic 
mode can now be obtained from (fT4|) . It will be con- 
venient to express it as D{q^) = —TT^^{q^), where 



d = 2 



d'^fc 



d'^fc 



fc2 



4mv 



iO+ 



(17) 



We have parametrized the two-fermion scattering length 
in the ground state by v and a mass scale m. Since the 
ground state is not vacuum but a completely filled valence 
band, resonant scattering may occur at ^ 0, which we 
shall determine later. We shall also exploit the freedom 
to define the mass scale m in the most convenient way. 

The expression (fTT]) is certainly valid in the saddle- 
point approximation employing the bare fermion prop- 
agators and bare vertices. However, its form does not 
qualitatively change as a result of quantum fluctuations 
in the insulating state. One can apply diagrammatic per- 
turbation theory to construct the exact bubble diagram 
and derive the exact boson Green's function. All excita- 
tions, bosonic and fermionic, are gapped, so the propaga- 
tors and vertices acquire only numerical renormalization. 
The exact fermion propagators contain dispersions (jl6p 
where the gaps Ac, Av and masses nic, have renor- 
malized values. While these quasiparticles and quasi- 
holes are mixtures of the microscopic particles and holes, 
their pairing rules are the same as in the saddle-point ap- 
proximation due to the Fermi-Dirac statistics. Similarly, 
since the renormalized interactions remain short-ranged, 
the only relevant consequence of vertex corrections is the 
renormalization of 7c and 7v. Therefore, the exact boson 
propagator is the same function of renormalized fermion 
parameters as ((T7)) . and we can use it to establish cer- 
tain relationships between exact bosonic and fermionic 
excitations in the insulating state. 



In two dimensions, integrating momenta up to a cut-off 
|A;| < A in (jl7p gives the unregularized expression: 



TOv7^ 

47r 



'•"c7c 

47r 



log 1 



A^/mc 



r2 - 2Ac - + ^o^ 



log 1 - 



A^/TOv 



-r! - 2Av - 4^^^ + ^0^ 



(18) 



In order to avoid an infra-red divergence, we must assume 
that fermionic excitations are gapped in the ground state 
when evaluating Iliog. Taking the particle and hole gaps 
to be ^c and 5^ in the ground state at = respectively, 
we find 

mz. mc7c% ^2Ac + 3^-^^-.0+^^ 



All Att 



Att 



(19) 



in the A — > oo limit. 

The boson propagator D{q^) = —T^{q^)~^ has branch- 
cut singularities at real frequencies $7 > Wc and D. < —ujv 
beyond the threshold for Cooper pair decay, where 

2 2 

Wc = 2Ac + , = 2Av + . (20) 

4mc Aniv 

The poles of D{qf^) sitting between these branch-cuts on 
the real axis yield the dispersion of charged bosonic ex- 
citations in the insulating state. In order to obtain them 
from the non- linear equation 7r(q^) =0, let us write 

mc7c — m + Sm , m^'-f^ = m — Sm , (21) 

which defines the masses m and 6m, and separate the 
terms proportional to m and Sm in (jl9p . We find that 
7r((7^) = is equivalent to 



(wc - Q) (cjv + ^) = ScS^e " ( ^ , ^ 



Therefore, the poles are obtained by solving: 



Wc - Wv , 
- ^ ± 



0' ■ 

(22) 



(23) 



First, we shall verify that the poles always exist by 
substituting (1^ into (ITOl) . After introducing 



h = 



2/3 



Wc + 



, bo 



AS^S^e-" 

(Wc + LO^Y 



(24) 



and some algebraic manipulation we obtain the condition 
equivalent to T^{q^i) = 0: 



(5v V^o 



ITVT^ 



1± 



(25) 
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FIG. 4. (color online) The existence of bosonic poles in two 
and three dimensions, (a) d = 2: A plot of (1 =p \/l — b) /(I ± 
y/T^) (solid) and (5c/5v) x (6/60)""'''^'" (dashed) from (gSl) 
for positive and negative 5m. bo — 0.7, 5m/m — ±0.1, 5c = 5^ 
in this example, (b) d = 3: A plot of Iie{n{q^)} — mi^/An 
from (|29p as a function of f2 G R. The shaded areas are 
sparmed by = Re{Q,}e''^ for 9 £ (-7r/2, 7r/2), the bottom 
boundary corresponds to 8 = 0. The two horizontal lines 
indicate —mvc/^n (upper blue) and — mj/v/47r (lower red). 
ojc/uJv — 0.5, {5m/m')^^'^ = —0.089, vertical axis is in arbi- 
trary units in this example. 



The FiglUJa) shows plots of the left-hand and right-hand 
sides of this expression. In general, there are two inter- 
sections between these curves for any &01 which can be 
regarded as a tunable parameter. One intersection cor- 
responds to the particle-like pole Vtc, and the other one 
to the hole-like pole Jlv. In the limit 60 the intersec- 
tions occur at small &, so that the right-hand side of 
behaves either as & or 1/5, depending on the sign choice, 
that is the particle/hole nature of the solution. Then, 
1 5m I < m guaranties that the curve on the left-hand side 
always intersects the functions b and 1/6 on the right- 
hand side when 60 —J' 0, regardless of the sign of 5m. In- 
creasing &o shifts the intersection points to larger values 
of b, and there is a "critical" value of 60 at which the two 
intersections coalesce and beyond which the poles disap- 
pear. However, the superfluid transition occurs before 
this "critical" 60 is reached. Namely, equating the pole 
frequency to zero at g = in (^5)) yields the condition for 
the superfluid transition: 



4AeAv = 



(26) 



which is equivalent to h — 4AcAv/(Ac -f Av)^ < 1. 
Therefore, the transition occurs before the value of b 
corresponding to either the particle-like or hole-like pole 
becomes 1, which by continuity is before 60 reaches 
the "critical" value. The limiting-case transition with 
& = 1 is special and involves condensation into relativis- 
tic bosonic modes. 

Knowing that bosonic charged excitations always have 
well-defined particle and hole-like poles in the insulating 



state, we can examine some of their important proper- 
ties. We can immediately see from that the bosonic 
spectrum is relativistic if fii = —Vl.2 at q — > 0, that is 
Ac = Av. Equal participation of both particles and 
holes is required for relativistic dynamics, and the re- 
sulting superfluid transition in the XY universality class 
occurs at the expected place. The particle-hole symme- 
try is generally lost at finite and especially large q, unless 
TOc = m-v. This is natural since the microscopic dynam- 
ics is not particle-hole symmetric, but this does not affect 
the XY universality of the transition, shaped by the long- 
wavelength modes. 

Finally, we observe that the boson gaps are always be- 
low the threshold for decay into fermion pairs. It imme- 
diately follows from ^ that 0.1 < 2Ac and jrJal < 2Av 
at q = for any /3. As a consequence, at any finite /3 
the superfluid transition f2 = occurs when Ac and Av 
are still both finite. Fcrmionic excitations are gapped at 
the transition and may be integrated out from the micro- 
scopic action to yield an effective bosonic action which 
describes the transition. The insulator adjacent to the 
superfluid is a bosonic Mott insulator. 



B. d = 3 

Following the same procedure as before and in Refi^, 
we find from (|17l) the regularized boson propagator 
D{q,,) = -7r{q,,)-^: 



niq,) = ^ -f ^y^ + -c(2Ac-f^)-^0^ 



TOv7v 1^ 



47r V 4 



+ m^{2A^ + n) + iO+ . (27) 



Again, there are branch-cut singularities for real frequen- 
cies n > ujc and < — Wv defined in ([^01) . Let us define 
mass scales m and Sm by: 

mij'^=m^+6m^ , mv^^ — m^—Sm^, (28) 

allowing Sm^^'^ to be either positive or negative, and then 
write 

, mV m I I — — I — — \ 

+ ^ (V"^('^c - f^) + ^Jm{uJs, + r?) j (29) 

v=v ^ ^ ^■\/m(a;c - il) - y^m{ujy + fi)^ . 

The solutions of 7r((7^) ~ define bosonic poles of either 
particle Re{f2} > or hole Re{57} < flavor. The con- 
ditions for the existence of the particle-like and hole-like 
poles at a wavevector q are v < vdq) and v < Vv{q) 
respectively, where 



VrXq)=- I 1 - ( — ) I ^H^c + c^v) (30) 



Vy{q)=- ( 1 + ( ) I \/m(Wc + 
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following from the requirement that the square roots in 
(|29p have positive real parts. An example is shown in 
FiglHb). The detuning parameter values v — i^dq) and 
v = Vv{q) correspond to resonant scattering between 
fermion particles and holes respectively, since they sepa- 
rate the regimes in which two-particle or two-hole bound 
states exist (BEC) or do not exist (BCS). When the poles 
exist they can be determined by solving 

(31) 

whose structure is the same as in (|23p . allowing us to 
immediately read out the condition (|26|) for the phase 
transition (set = and g = 0), 

4AcAv = (^Ae -f Av - ^) . (32) 

The transition becomes pair-breaking if it occurs at a 
vanishing fermion gap Ac = or Av — 0. The former 
requires v = —y/2rnK^ or equivalently v — J^c(O), and 
the latter requires v = —y^2mAc oi v — i'v(O). 

In the BCS regimes there are no bosonic poles, so that 
the only path to superfluidity involves the standard BCS 
instability when a fermion gap is closed to produce a 
finite quasiparticle density in the ground state. In the 
BEC regimes, real poles exist in the insulating state (the 
occurrence of complex poles in ((3T|) is preempted by a su- 
perfluid transition) . Both the quasiparticle and quasihole 
gaps, Ac and Av respectively are finite at the transition 
in the BEC regimes, so that this transition is in a bosonic 
universality class. 



V. DISCUSSION AND CONCLUSIONS 

Attractive interactions between fermionic quasiparti- 
cles which remain gapped across the superconducting 
transition give rise to a protected bosonic Mott insula- 
tor of Cooper pairs. The superconducting transition at 
a fixed density is associated with phase fluctuations of 
the order parameter, and the insulating state adjacent 
to the superconductor is characterized by gapped low- 
energy collective modes with bosonic statistic. In two 
dimensions this Cooper pair insulator completely sur- 
rounds the superconducting phase in the phase diagram, 
making a direct transition between the superconductor 
and an unpaired band insulator impossible. Above two 
dimensions, the bosonic Mott physics belongs to the BEC 
regime of sufficiently strong coupling, while the transition 
is pair-breaking in the weak coupling limit. Although suf- 
ficiently strong pairing and perturbations can stabilize in- 
sulating symmetry broken phases of Cooper pairs, there 
is no a priori reason for this Mott insulator to be ther- 
modynamically distinguished from the band insulator. A 
sharp distinction nevertheless exists in non-equilibrium 



settings, and can be characterized as a non-equilibrium 
phase transition. 

The most abstract and unbiased way of identifying 
transitions must be representation independent. Such in- 
formation is contained only in the energy spectrum of the 
Hamiltonian, or density matrix as a generalization which 
can describe coupling to the environment. Non- analytic 
changes of the total density of states in the ground state 
manifold are associated with thermodynamic quantum 
phase transitions. Analogous non-analytic changes in the 
excited state spectrum are possible, and can show up as 
sharp phase-transition-like phenomena in systems driven 
out of equilibrium. We demonstrated that the presence 
of Cooper pairs in the bosonic Mott insulator can at least 
be identified with such a non-equilibrium transition, and 
physically associated with the ability to drive Bose con- 
densation among excited states into a Cooper pair laser, 
the state analogous to the many-body state of photons 
produced by a laser. 

We considered a band insulator subjected to pairing as 
a model system. The simplest realization of such a sys- 
tem is found in trapped neutral ultra-cold gases of alkali 
atoms placed in an optical lattice. The density of atoms 
can be chosen to correspond to two atoms per lattice 
site in the central portion of the trap, while the strength 
of attractive interactions among them is routinely con- 
trolled by the Feshbach resonance. A superfluid transi- 
tion from a thermally excited band insulator has been 
already experimentally studied in this kind of a system 
in the vicinity of the BCS-BEC crossover^. Detecting 
a Bose-insulator behavior may be more difficult, as it 
requires very low temperatures and a way to sharply dis- 
tinguish the presence of Cooper pairs (locally stable pairs 
have been observe d^^i^"^ ). Perhaps the excited-state con- 
densation might be possible to drive by an external low- 
frequency radiation, and then detect in the momentum 
distribution of particles in a time-of-fiight measurement. 

In condensed matter systems, a two-dimensional band 
insulator could be rendered unstable to Cooper pairing 
by the proximity effect to a superconducting layer (in a 
manner similar to Ref.~-). The challenge here is to find a 
material whose bandgap is comparable to the typical crit- 
ical temperature of the superconducting layer. Bismuth 
antimony alloys seem to be promising candidates, since 
their bandgap can be tuned by the antimony concentra- 
tion between zero and at least 25 meV=. Another solu- 
tion might be to exploit a charge-density- wave or an engi- 
neered superlattice to create a small bandgap. Detecting 
bosonic excitations in such insulators also seems difficult, 
but indirect indications might be observable in transport 
measurements, such as Nernst effect in the vicinity of 
the superconducting transition. The Cooper pair insu- 
lator physics might also be relevant for the amorphous 
oxide films. 

Perhaps the most promising way to test the Cooper 
pair insulator hypothesis is by Monte-Carlo simulations 
in negative-fj Hubbard models^. External driving mech- 
anism can be simulated, and it should be possible to 



13 



observe the formation of a driven condensate, and the 
corresponding non-equihbrium transition to a true band 
insulator. 



A. A qualitative picture of cuprates 

The apphcabihty of the picture presented in this paper 
to cuprates depends on two things: (i) whether fermionic 
excitations acquire at least a partial gap independent of 
pairing®^, for example due to spin correlations or fluc- 
tuations in the particle-hole channel, and (ii) whether 
we can understand the two-dimensional c?-wave super- 
conductivity in the same spirit as a BEC regime despite 
some obvious differences, like the absence of electron pair 
bound states. 

Phenomenologically, an independent antinodal gap in- 
deed exists and is reflected by the T* crossover tempera- 
ture in underdoped cuprates which grows the closer one 
gets to the antiferromagnetic (AF) phase. In the context 
oi t — J models one could naively relate the fermion gap 
to the exchange energy cost ~ J{S) for removing a spin 
from the local AF environment, where (S) is the effec- 
tive local magnetic moment per site, reduced by quan- 
tum fluctuations due to doping. There are a number of 
sophisticated theories dealing with the microscopic ori- 
gin of such a gap^"^. For our purposes, it is sufficient 
to assume that the gap originates from spin dynamics 
and that it decreases with doping. This would ensure an 
antinodal gap dominated by the nearest-neighbor spin 
exchange. 

Next, we recall an old idea that charge organized in the 
form of Cooper pairs could propagate through a short- 
range correlated AF environment without frustration. A 
familiar caricature is that of a hole moving by nearest- 
neighbor hopping and leaving a trail of frustrated spins 
along its path; a second hole which follows right be- 
hind would repair this damage and thus the two are " at- 
tracted" to each other by the intervention of the AF en- 
vironment in which both move. Some variation of this 
mechanism - whose precise form is unknown at present 
time - then suffices to produce an effective attractive 
short-range interaction between holes, which unavoidably 
has anisotropic properties, being the strongest along the 
lattice bond directions. Certainly, this interaction has 
to fight against "quantum entropy" and, at a given low 
doping, need not be strong enough to actually cause pair- 
ing, let alone superconductivity. However, its anisotropy 
ensures that antinodal holes are most affected. 

Therefore, we arrive at a situation which is similar to 
the one analyzed in this paper. Fermionic charged exci- 
tations near antinodal regions of the Brillouin zone are 
gapped, but experience short-range (anisotropic) pairing 
interactions. If the concentration of holes is increased, 
the gap decreases, while the effect of interactions be- 
comes stronger because the average separation between 
holes shrinks. These trends persist until the character- 
istic interaction scale becomes comparable with the gap 



scale, and then the conditions for superconducting in- 
stability are created. The analysis in this paper shows 
that bosonic excitations with charge 2e and no spin nec- 
essarily play a crucial role in these circumstances. The 
short-range "antinodal" Cooper pairs have a certain ro- 
bustness across the superconducting transition and pro- 
vide a backbone for a correlated "normal" state. 

However, the d-wave order parameter symmetry re- 
quires that nodal quasiparticles remain gapless across the 
transition, which makes cuprates fundamentally differ- 
ent from the models studied in this paper. A conven- 
tional BEC regime is out of question since it requires a 
full quasiparticle gap and a Cooper pair bound state. It 
is useful to think of the difference between the two in 
the real space. Deep in the conventional BEC regime, 
tightly bound Cooper pairs have a small typical size, 
smaller than the average separation between particles, 
and the pairs much larger than this size are absent for 
all practical purposes. In the nodal d-wave case there is 
a power-law distribution of the pair sizes and thus there 
are arbitrarily large Cooper pairs, testifying to the pres- 
ence of nodes. Nevertheless, we expect that in cuprates 
most of the Cooper pairs are still rather small, their size 
basically given by the Cu-O-Cu bond length. It is these 
small, stick-like Cooper pairs that, within our picture of 
the pairing gap, are responsible for the antinodal gap. 
These bond Cooper pairs can actually be represented 
by an effective bosonic theory of the type discussed in 
this paperiii While such purely bosonic theory is in- 
complete when applied to cuprates and must be supple- 
mented by the nodal fermionic quasiparticles, a subset of 
its properties in the charge sector is in fact shared by the 
tightly bound BEC-limit pairs. For example, there will 
be preference for certain "magic dopings" where Cooper 
pairs can neatly fit into the underlying lattice, although 
the resulting "pair density wave" patterns will in gen- 
eral be rather different - the reason is that, even at this 
level, there is a difference from s-wave description since 
the stick-like nature of d-wave Cooper pairs in cuprates 
makes them unusually prone to formation of various quasi 
one-dimensional and nematic patterns ordinarily absent 
from the conventional BEC limit Very importantly, the 
decay channel for Cooper pairs in cuprates is anoma- 
lously slow due to the linearly vanishing density of quasi- 
particle states at the gapless points in the Brillouin zone. 
The resulting longevity of Cooper pairs is what likely 
protects the integrity of bosonic fluctuations in the pseu- 
dogap state, both in terms of pairbreaking and vortex 
dynamics. "^^ 

Thus, we hypothesize that charged quasiparticles in 
cuprates live in an unconventional BEC-like regime due 
to both quasi-two-dimensionality and relatively large 
strength of interactions. This BEC regime is peculiar 
since it does not involve low-energy bound-state pairs as 
a result of the d-wave pairing symmetry. Consequently, 
it lacks a true, sharp charge gap and it clearly has gap- 
less nodal spin excitations. While some of the features 
of this picture can be gleamed from the d-wave dual 
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approach,-S4 many details needed to fully verify its va- 
lidity are yet to be worked out. For example, the emer- 
gence of this peculiar BEC regime from within a renor- 
malization group analysis^^ might require non-analytic 
features of the effective interactions between quasiparti- 
cles at zero momentum transfer. In any event, this may 
be the route to ultimately understanding a variety of un- 
conventional properties of cuprates as manifestations of 
quantum vortex dynamics, including the Nernst effect 
and density-wave patterns observed by Scanning Tun- 
neling Microscopyiiii^i^i^. 

At very low dopings, close to the AF phase, the domi- 
nant low-energy excitations are bosons which carry spin 
5=1, but no charge. Therefore, a certain crossover 
between spin-dominated and charge-dominated regime is 
expected in the pseudogap state at zero temperature. It 
is natural to interpret the region in the cuprate phase 
diagram where the Nernst signal is large as a charge- 
dominated regime, below a dome-shaped onset tempera- 
ture for pairing which encloses the superconducting dome 
in the phase diagram of cuprates. 
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Appendix A: Effective theories of pairing in band 
insulators 



Here we derive effective low-energy bosonic lattice the- 
ories from the microscopic continuum model ([1]) in the 
regime where the lowest energy excitations with bosonic 
statistics lie well below twice the finite gap for fermionic 
excitations. We can introduce a complex Hubbard- 
Stratonovich field <I>(ri,r2) to decouple the interaction 
term in ([l}: 

-$(ri,r2)t$(ri,r2) 



J -]d 



U{\ri -r2\ 
-F$(ri,r2) Vt('^i)V'l(r-2) + /1-C. 



(Al) 



In order to rewrite the model in a tight-binding form, we 
represent the field operators as expansions weighted by 



Wannier functions: 



(A2) 



nM 



ni Ri ,n2-R2 

rii,Ri 712, R2 

Here, R are the discrete coordinates of lattice sites, 
n band-index or orbital labels, fanR fermionic lattice 
fields and 6niKi,n2R2 bosonic lattice fields. The bosonic 
Wannier functions W can be chosen at will (subject to 
constraints discussed below), while the best choice for 
the fermionic Wannier functions W is dictated by the 
bare fermion dispersion in the lattice potential. If Tpnk 
are Bloch wavefunctions in the potential V{r), then the 
choice 
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(2^ 



-ikR 



(A3) 



allows us to express the fermionic kinetic part of the ac- 
tion as 



dr 



fanR 

.nR 



d_ 



LnR m 



EE 
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n R1R2 



where the hopping couplings are functions of the bare 
energies e„fc: 



tnAR 



The quadratic bosonic term becomes: 
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where 

KniRi,n',R\,n2R2.nLR'2 ~ / '^'''^l '^'^'"2 TTT] TT X 

W^„*,R,(n)W^„;K;(r-i)W?4^,(r2)W^„,Hi(r2) . (A7) 

The localized nature of Wannier functions will generally 
result with significant values for K only when Ri R!^ 
and R2 ~ i?2. On the other hand, the interaction po- 
tential U{\ri — r2\) controls the dependence of K on 
AR = Ri- R2 and AiJ' = R!^ ~ R'^, causing K to 
become larger and larger as Ail and Ail' grow. Since 
K acts like a mass term for the b fields, fiuctuations of 
long-range singlet bonds will be suppressed. The require- 
ment that K be convergent puts certain restrictions on 
the choice of W. Namely, any short-range interaction po- 
tential ?7(|ri — r2|) will generate diverging factors in the 
limit Iri — r2| — > oo which must be controlled and tamed 
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by sufficiently fast decaying Wannier functions. With the 
exception of the contact potential U (x 5{ri — which 
we will separately discuss later, we shall assume for sim- 
plicity that U{\ri — r2\) is not too short ranged, so that 
we can use the same Wannier functions for the boson and 
fermion fields, WnR.{r) — WnR{r). Then, the interaction 
term takes a simple form: 



Sbf — / ^ ^ (^niRi,n2fi2-^t«lfil'^in2R2 + 



h.c. 



niRi n2R[ 



(A8) 

The full lattice action S = S f + Si, + S^f is an exact 
rewriting of the starting model We can introduce 
a short-hand notation for (n, R) by defining a "hyper- 
lattice" whose sites i,i ■ ■ ■ carry both the regular lattice 
coordinates R and band/orbital coordinates n: 



S= dT 



ijkl ij 



(A9) 



Since we assume that the chemical potential fi lies in a 
bandgap, the bare fermionic fields are fully gapped and 
safe to integrate out. We are interested in the low-energy 
dynamics of singlet bonds or Cooper pairs bij , so we also 
integrate out all high energy bij fields obtained for too 
large spatial separations between sites i and j. The ef- 
fective theory containing only the degrees of freedom be- 
low the fermion bandgap scale is a purely bosonic lattice 
model in which the b fields acquire effective dynamics. 
We cannot derive the exact form of the effective action, 
but schematically: 



ijkl 



d_ 



^kl 



(AlO) 



ijkl i'j'k'V 



Since only the short-range bonds bij remain, the num- 
ber of bosonic degrees of freedom scales as the number 
of lattice sites, with an addition of a finite number of 
low-energy "orbital" states per spatial lattice site. With- 
out knowing the couplings, this action is very general. 
We can be only certain that it captures a superfluid- 
insulator transition, because the "two-channel" model 
(lAip is guarantied to prefer the condensation of $ at any 
U^^ when /i approaches close enough to a band-edge. 
Since "orbital" states are separated by finite gaps, only 
one bosonic mode becomes soft at the superfluid tran- 
sition, implying that a different bosonic effective model, 
with one degree of freedom per site, is also sufficient to 



describe the universal properties of this transition. Be- 
low we outline the derivation of such an effective model, 
which discards all microscopic details of the interaction 
potential J7(|Ar|) and replaces it by a contact interac- 
tion. However, keeping these microscopic details in the 
action may be needed for the description of other transi- 
tions between the insulating states. For example, should 
in some circumstances the singlet bonds become strongly 
repulsive, yet close-packed, the effective theory of this 
kind would be able to describe a variety of Mott insulat- 
ing states, including valence-bond states. 

If one is mainly interested in universal properties, one 
identifies a renormalization group fixed point and con- 
siders only the scale dependence of relevant operators. 
The microscopic details of the interaction potential U 
are found to be irrelevant at least at the Gaussian fixed 
point, and at unitarity [/ — ^ cx) in (jAl[) in dimensions 
d > 2. The only relevant piece is the contact interaction 
C/(|Ar|) — ^ US{Ar), for which the above bond effective 
theory breaks down. Instead, one may decouple the inter- 
action by a Hubbard- Stratonovich field $(t'), and define 
on-site bosonic fields b by: 



$(r) 



E WnR{r)br 

n,R 



R 



(All) 



The full lattice theory takes the form: 



i ij k 

and the effective pure bosonic theory is: 



, (A12) 
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Appendix B: The range of effective interactions 
between bosonic excitations in the Mott insulator 



Here we outline the calculation of the diagram in 
FigHJb) which describes the scattering processes of 
charged bosonic excitations. We assume that the ground 
state is a bosonic Mott insulator so that coherent bosonic 
excitations exist and all excitations are gapped. We 
are interested in low energy bosons which cannot decay 
into fermions during the course of a collision. The bo- 
son and fermion gaps are Ab and At respectively, where 
Ab < 2Af. 

Following Feynman rules and writing = {uj,k), 
Qif^ = (r2,;,qj), = we start from the generic 

scattering diagram in Figlljb): 
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Assuming that all fermionic excitations are particle-like, 
with energy 



f.2 

Af + — > 
2m 
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each denominator in the second line contains the same 
implicit infinitesimal tcTin ?'0+. The internal frequency 
integral extended to a loop over the upper complex half- 
plane picks the poles of the last two Green's functions. 
We have also introduced 



ai = fii — £_ 



Q!2 = ^2 - e-fe+q^+f + 2 ■ 



for notational brevity, and we will use S-k 
convenient. 
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Next, we concentrate on zero energy f2 = and fi- 
nite momentum p = \p\ transfer. This will allow 
extracting the range of static (non-retarded) effective in- 
teractions. Also, let us consider only the lowest energy 
incoming particles whose momenta arc zero, Qi = Q2 ^ 0. 
In order to keep these incoming particles on the mass 
shell we must also choose fii = O2 = Ab. Then: 
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where 5 = 2Af — Ab. It is already clear that 4m(5 is 
the only externally set scale which appears alongside . 
Since the Fourier transform Uh (p) determines the spatial 
dependence of the effective interaction potential, we see 
that r ~ {4:m5)~^^^ is a length-scale characterizing this 
interaction. As long as Ab < 2Af, r is finite and the 
effective interaction is short-ranged. 



Further progress can be made in d 
dimensions. Changing variables to k = 



= 2 and d = 3 
k'^/m we obtain 



in d = 2: 
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This expression is finite for any p as long as (5 > 0. In 
(i = 3 we can write x = cosd and immediately integrate 
out the polar angle: 
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where in the last step we have taken the limit <C 4m(5 
before integrating out k. Interpreting these expressions 
as Fourier transforms of the effective interaction potential 
we clearly find that in both two and three dimensions the 
effective interaction has a finite range r ^ (4m(5)"^/^. 

If incoming particles had finite momenta and ener- 
gies above the boson gap, no qualitative modifications 
of these conclusions would be encountered other than 
i5 — >■ 2 Ay — e, where e can be the energy of either in- 
coming boson. This could be interpreted as an energy- 
dependent effective interaction range, but it also reflects 
the fact that too energetic particles break into fermion 
pairs during collisions, so that scattering cannot always 
be described merely by an interaction potential. 
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